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Abstract 

The effect of smooth shape changes of metahic nanoparticles on localized surface plasmon res- 
onances is assessed with a boundary integral equation method. The boundary integral equation 
method allows compact expressions of nanoparticle polarizability which is expressed as an eigen- 
mode sum of terms that depends on the eigenvalues and eigenfunctions of the integral operator 
associated to the boundary integral equation method. Shape variations change not only the eigen- 
values but also their coupling weights to the electromagnetic field. Thus rather small changes 
in the shape may induce large variations of the coupling weights. It has been found that shape 
changes that bring volume variations greater than 12% induce structural changes in the extinction 
spectrum of metallic nanoparticles. Also the largest variations in eigenvalues and their coupling 
weights are encountered by shape changes along the smallest cross-sections of nanoparticles. These 
results are useful as guiding rules in the process of designing plasmonic nanostrucrures. 

PACS numbers: 41.20.Cv, 71.45. Gm, 73.20.Mf 
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I. INTRODUCTION 



Localized surface plasmon resonances (LSPRs) have their origin in the interaction of 
metalhc nanoparticles (NPs) with hght. The advancements made in the last decade have 
enabled the use of LSPRs for light manipulation at the nanoscale^. The applications of 
LSPRs include enhanced sensing and biosensors^''^, cancer imaging and therapy^, plasmonic 
lasers^ and spasers^, enhanced nonlinearities, etc. Metallic NPs can be made by chemical 
synthesis or by lithographic techniques. Thus a large variety of shapes have been obtained 
by chemical synthesis. Beside the spheres, the shapes include nanorods^, cubes^'^, triangular 
prismsi^, tetrahedra^, and hexagonal prisms^.. In contrast, in "top-down" or lithographic 
techniques the shape of the metallic NPs are "flatter" such that structures like disks^i, 
dimers of diska^^ or bowtie structures^ have been successfully fabricated. 

Optical properties associated with LSPRs are determined by the shape, size, structure, 
and local dielectric environment of the NPsi^^i^. Along with the electric field enhancement 
created around NPs, the effect of the dielectric surrounding the metallic NPs is widely 
used in sensing by monitoring the shift of LSPR absorption peaks with respect to the lo- 
cal dielectric changes. The absorption peak redshifts as the embedding refractive index is 
increased^ii^.. On the other hand, in the quasistatic approximation, when the particle size 
is much smaller than the wavelength, the size of particle doesn't play any role^^. Thus, 
in the quasistatic approximation, a metallic nanosphere exhibits just a single LSPR which 
is dipolar, irrespective of the size. However, as the radius is increased and the quasistatic 
approximation is no longer valid, the LSPR of the nanosphere shifts toward infrared and 
higher multipolar resonances emerge in the spectrum^^^ii^. Moreover, by elongating or flat- 
tening the spherical NPs, one obtains spheroids which have two LSPRs corresponding to 
longitudinal and transverse polarization of light. For spheroids the plasmon resonance re- 
sponse depends solely on the aspect ratio, which is the ratio between the lengths of rotation 
axis and thickness of the particle^*^. The prolate spheroids have the aspect ratio greater 
than 1, such that the longitudinal mode shifts to longer wavelength, whereas the transverse 
mode shifts to shorter wavelength. On the contrary, in oblate spheroids the longitudinal 
mode shifts to shorter wavelength, whereas the transverse mode shifts to longer wavelength. 
Thus by simply adjusting the aspect ratio, the LSPRs can be tuned at will over a broad 
range of wavelengths. 
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The discrete-dipole approximation (DDA)^ the finite- difference time domain (FDTD) 
scheme^, and the boundary element method (BEM) in a full electromagnetic calculation^^ 
have been successfully used in the calculation of the optical response of arbitrarily shaped 
NPs. These methods integrate the full Maxwell's equations but they are numerically ex- 
tensive. Thus they cannot be directly used in the problem of designing plasmonic nanos- 
tructures. In addition, the above methods offer little inside about the formation, nature, 
and the behaviour of the LSPRs. As a more physical approach that works very well in the 
quasi-static limit, the hybridization model^i has been proposed to solve some of the above 
issues regarding DDA, FDTD, and BEM methods. On the other hand, also in the quasi- 
static limit, the boundary integral equation (BIE) methocy^i^I enables a direct relationship 
between the LSPRs and physical parameters of NPs like the shape (geometry) or the com- 
plex dielectric functions by relating the LSPRs to the eigenvalues and eigenf unctions of the 
operator associated with BIE. 

In the current work, the LSPR spectral modifications made through small but smooth 
shape variations from spheroids are studied by the BIE method of 16. Smooth deformations 
from spherical shape, prolate, and oblate spheroids are considered by keeping the same as- 
pect ratio. The parameter that uniquely describes the shape variations can be related to the 
relative volume variation from spherical shape and from prolate and oblate spheroid, respec- 
tively. The infiuence of small shape changes on LSPRs has been treated in several works 
with emphasis on NP roughness^Si^— or smooth shape variational. All these papers monitor 
the spectral shift and, eventually, the inhomogeneous broadening of the main resonance due 
to shape variation. Recently a perturbative method has been developed in order to calculate 
the eigenvalue changes of LSPRs at small shape perturbations^. Despite many advantages 
like the use for designing plasmonic nanostructures with predetermined propertiesSE^ the 
method does not provide directly the weights of the LSPR eigenmodes and their modifica- 
tions when small shape changes occur. On the other hand, in addition to the fact that it 
can be adapted to the recipes of^*^, the BIE method provides both the LSPR eigenvalues 
and their weights^^. The present work shows that not only the eigenvalues change with the 
smooth modifications of the shape but also their weights change, sometimes in a drastic 
manner. There are also other goals of the present study. One other goal is related to the 
nanoparticle design and fabrication, which must be robust against the variation of nanopar- 
ticle shape. Thus during the fabrication processes like chemical synthesis or "top-down" 



3 



approaches based on e-beam lithography variations of the process parameters are encoun- 
tered. An example is the lift-off step in the top-down approach, where some precautions have 
to be made in order to fabricate systems with small features like dimers^S or oligomers^i. 
The present paper shows numerically that the eigenvalues and their weights have the largest 
variations when the applied field is along the smallest cross-sections. Hence additional care 
has to be considered in order to fabricate successful metallic NPs. Another reason comes 
from the simulations of experimental data, where one must find the most realistic model 
that fits the experimental setup. Thus in the simulation process some spectral features may 
be attributed only to the retardation and the shape variations are not considered at all. 
Generally, retardation moves the electrostatic resonances toward infrared^ii^. At the same 
time, some eigenmodes that are dark in the quasi-static limit become visible when the re- 
tardation is considered. In the present study it is shown that redshifting and the emergence 
of other higher order LSPRs may be obtained also by small shape changes. 

The paper has the following structure. The next section presents the BIE method and 
its accuracy. Then, the LSPR modifications due to smooth changes from various spheroidal 
shapes are analyzed in the section that follows the next section. The last section is dedicated 
to conclusions. 

II. THE METHOD AND ITS ACCURACY 

In the electrostatic (quasistatic) approximation the optical behavior of the metallic 
NPs is described by the Laplace equation, whose solution may be obtained by the BIE 
method^'^^^. The metallic NP delimited by surface S is assumed to have a complex 
permittivity ej and is immersed in a medium of complex permittivity eg. The incident 
electromagnetic field is represented as a uniform electric field Eq. The electric potential $ 
associated with the total field E obeys the Laplace equation A$ (x) = 0; x G with 
boundary conditions eoff |_,_ = x G S and — V$ (x) — Eq, |x| — )• oo. The total 

electric field is E = — V$ (x), n is the normal to the surface S, and 3?^ is the 3D Euclidian 
space. The electric potential can be expressed in terms of the single-layer potential the 
surface charge fiEo as 
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<l>(x) = -x-Eo + -^ / (1) 

An J |x-y| 

x.yGS 

The charge density fiEo obeys the following integral equation 

^/i£;Jx)-MH =n-Eo, (2) 
where n • Eg is the dot product of vectors in 3D, M is defined on surface S as 

x,ygS 

and A = (e, — Cq)/ (cj + eg) is a dielectric factor that depends on eo and ei. The operator M 
and its adjoint 



x,yeE 

have the same discrete and real spectrum that is bounded between 1/2 and -1/2. In addition 
the number 1/2 is an eigenvalue irrespective of the NP shape. The charge density is 
provided in terms of the eigenvalues and eigenfunctions of M and as 

/^i?o = XlxZ — \uk) {vk\ n ■ Eo) . (5) 

k 2A 

where \uk) , |ffc) are the eigenfunctions of M and M^, respectively, Xfc is the k*'' eigenvalue of 
M and M''", and n ■ Eq) is the scalar product of two square- integrable functions defined 
on E. The specific polarizabiliy of the NP that is the dipole moment generated by fiEo 
divided by the NP volume V has the form of a sum over all eigenvalues of M and 

k 2A Afc 

In Eq. (E]) Pfc = (x ■ N| Uk) {vk\ n ■ N) /V is the coupling weight of the k^^ eigenmode to the 
electromagnetic field and N is the unit vector of the applied field given by Eq = EqN. The 
coupling weight pk is defined in terms of the eigenvectors of the operator M and its adjoint 
M^and it depends only on the shape of the NP, while the specific polarizability a depends also 
on the dielectric permittivities^. The eigenmodes with a non-zero pk couple with the light 
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and they are called bright modes. On the other hand, many other eigenmodes do not couple 
with light and have vanishing pkS, thus they are called dark modes. Equation ([6]) shows 
an eigenmode decomposition and allows an analytic expression for specific polarizability a 
when the dielectric permittivity of the metallic NP is described by a functional form like 
the Drude model e (w) = £m — ^^/(u; (w + 27)). Here Em is the interband contribution to 
the permittivity including £00, i^p is the plasma frequency, 7 is the dumping constant, and 
i = Thus for an embedding medium with a real and constant dielectric permittivity 

Ed, the specific polarizability a in the Drude-like approximation takes the following form 
(Sandu et al. 2011) 

, , Wk{Em — ^d) Wk Ed ^pk /«x 

a[uj) = y ^ r^, (7) 

where 



'^pk - — i«J 

^eff.k 

is the square of the localized plasmon resonance frequency in the limit of negligible 7 and 
^eff.k = (1/2 + Xk) £d + (1/2 — Xk) £m IS an effective permittivity. The term (1/2 — Xk) is 
the depolarization factor (DF) of the eigenmodei^»2S. The specific polarizability a plays its 
role into the LSPR spectrum by its imaginary part that is proportional to the extinction 
spectrum of the incident lightr^. 

To evaluate the polarizability a one needs to calculate the eigenvalues and the eigenfunc- 
tions of M and using a basis set of functions that are defined on surface S. If the surface 
S can be parameterized by {x = g{z) cos (f,y = g{z) simp, z}, where g{z) is a smooth but 
otherwise arbitrary function, the variables that determine the surface are the coordinate z 
and the angle ip, respectively. This parameterization allows the surface S to be smoothly 
mapped onto unit sphere. Here g (z) is sufficiently smooth defined on z G [zmin, Zmax], such 
that g (zmax) = g {zmin) = 0, and g' (zmax) = g' (zmin) " ^ oo. By a linear transformation of 
the form z = ^ {{zmax — Zmin) zi + Zmax + -^mm), the ucw Variable zi is restricted to [—1, 1], 
such that the one-to-one mapping onto a sphere becomes now apparent. Thus, through 
the above mapping, the basis of spherical harmonics Yjm (6', p) defined on the unit sphere 
generates the basis of functions on S. All the other details about numerical procedure can 



be found in 
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FIG. 1: The comparison between the numerical results (the black dots) and the analytical results^ 
(the red full curve) of the DF of spheroids with various aspect ratios. 

The numerical accuracy of the method will be verified below for spheroids of various aspect 
ratio values. For ellipsoids and, in particular, for spheroids, there are three equally weighted 
bright modes determined by the electric field polarizations that are parallel to each axis. 
These bright modes are, in fact, the dipole eigenmodes, of which two of them are degenerate 
for spheroids. In addition to that the DFs of spheroids can be expressed analytically as 
functions of the aspect ratio defined by the term a^/a^ (spheroids being considered to be 
the ellipsoids of semiaxes a^, ay, a^, such that = dyY^- Thus, for spheroids, the following 
relations must be fulfilled: Px = Py = Pz = ^, Xx = Xy^ and theXx + Xy + Xz = 1/2 (Bohren 
and Huffman 1998; Sandu et al. 2010). In Fig. [T]the calculated DFs are compared to the 
analytical expressions provided in IS. A relative accuracy of at least 10~^ is obtained by 
using a basis of only 25 functions and a numerical quadrature of 96 Gauss points^^. Axially 
symmetric NPs permit much faster and more stable and accurate solutions due to the fact 
that numerical integration is performed only in 2D^~— . In fact, axial symmetry allowed us 
to obtain stable bright eigenmodes whose eigenvalues are close to 1/2, which is the largest 
eigenvalue of the system and is also dark— 
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III. THE EFFECT OF SHAPE VARIATION 



A. The shape 

Since the spheroids have an analytical form for their DFs, the variation of the LSPR 
response with respect to the aspect ratio change can be easily assessed. In the present work 
smooth shape changes that keep the same aspect ratio will be considered. Thus the NP 
shapes that will be taken into account below have the form 



giz) = Ciz,b)^l-[^y. (9) 

The parameters a and b, and the smooth function C {z,b) establish the shape of an 
individual particle. Without C {z,b), Eq. ([9]) describes a spheroid with an aspect ratio 
given by a and Zmax = —Zmin = Since the results of the quasistatic approximation are 
scale-invariant, just one parameter suffices to define a spheroid, which is usually described 
by two parameters (i. e., its distinct semiaxes). The function C {z,b) quantifies the smooth 
deviation from the spheroidal shape and is chosen to be peaked in the vicinity of 6 = 0, such 
that the spheroid is determined by 6 = 0. Thus the following form of C {z, b) is adopted 

Equations ([9]) and (ITO!) cover a large class of shapes including shapes similar to nanodisks 
and nanorods. For example, a nanorod with an aspect ratio of 4 can be generated by taking 
b=0.04 and a=4 and a nanodisk with an aspect ratio of 1/2 can be made by having b=2 
and a=0.5. Moreover, instead of parameter 6, one can choose a more intuitive parameter, 
the relative volume variation from a spheroid, as a shape parameter. Thus relative volume 
variations ISV jV are negative/positive for negative /positive b. 

In the calculations of the NP polarizability only the bright modes contribute to LSPR 
spectrum. Smooth deformations from spheroidal shape given by Eq. ([9]) will determine also 
higher order eigenmodes to become bright. For some of these bright modes, however, their 
weights may not be large enough to be distinguishable from the background or from other 
modes due to plasmon ovelapping. If 7is much less than the distance between two consecutive 
plasmon frequencies Wpfc? then the overlapping is not encountered and the following criterion 
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may be used to determine if a bright mode shows up into the LSPR spectrum. Eq. ( fTTl) can 
be used as a guiding criterion for the relevance of an eigenmode to the LSPR spectrum. It 
says that the height of the LSPR should be larger than a baseline of magnitude one. Thus 
for a ratio Copk/l ~ 10; which is perfectly normal, Ed/ ~ 1 (also within a typical range), and 
1/(1/2 — Xk) > 1, the weights fulfilling the condition pk > 0.1 would allow (fTTj) to be true 
and the bright eigenmodes to be observed in the spectrum. 

B. Changes due to the variations from spherical shape 

Throughout this work I consider small shape variations of gold NPs immersed in water 
('^i«ater=l-33). The Drude parameters of gold are those that are usuahy used in the literature. 
Thus, in energy, the plasma frequency cOp is about 9 eV and 7 is about 0.1 eV— . Also the 
interband Em is adjusted to have the LSPR for a nanosphere at 530 nm, i.e., Em = 10.2. 
Fig. [2] illustrates the variations of the two most relevant eigenvalues (or equivalently the 
associated DFs) and their weights with respect to shape variation from a nanosphere and field 
polarization. An eigenmode is considered relevant if its weight may exceed 1%. The relative 
volume AV /V varies with respect to spherical shape from -14% to 10%. Fig. [3] presents the 
variation of LSPR spectra with respect to both shape variation and field polarization. The 
shape variations are hardly detected by simple vizualization as one can see from the inset 
of Figinb, which shows the variations in cross-section for volume variations up to 14%. 

For parallel polarization, the first relevant eigenmode represents the dipole response and 
has a monotonic ( i. e. increasing) behavior of its DF with respect to volume variation. 
Therefore, according to Eq. 7 a smaller DF means a redshifted LSPR as it is observed from 
Fig[2K and [3^ at negative volume variations. The NP is invariant with respect to coordinate 
change z -H- —z, thus the quadrupole eigenmode is not allowed to be bright, therefore the 
second most relevant eigenmode is the octopolar eigenmode. In contrast to the dipolar one, 
the octopolar eigenmode exhibits a non-monotonic variation of the DF with a minimum 
reached when the particle is spherical {AV/V = 0%). No major variation of its weight (i. 
e. no more than 10%), however, occurs with no noticeable change in the LSPR spectrum. 
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FIG. 2: Changes of the two most relevant DFs (left side and continuous lines) and of their weights 
(right side and dotted lines) with respect to shape variation from sphere. The shape variation is 
accounted by the volume variation AV/V. The first (dipole)/the second (octopole) eigenmode is 
denoted by the square/star symbol. The field polarization is (a) parallel or (b) perpendicular to 
the rotation axis. 

Thus, from Fig [3^1 and according to the first sum rule X^Pfc = 1 (Sandu et al. 2010), the 

k 

weight variation of the dipole eigenvalue is also less than 10%. 

There is a different picture for transverse polarization (Fig. [2)d and FiglSb). The DF of 
the dipole eigenvalue is non-monotonic with a maximum around AV/V = —10 % and its 
weight varies by -20% for a volume variation of -13.6%. Moreover, for positive AV/V, the 
LSPR of the dipole eigenmode redshifts as a result of the DF drop. The large reduction of 
the dipolar weight is reflected in large increase of the octopolar weight, which can be also 
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FIG. 3: The LSPR spectra for three different AV/V: -0%, -7%, and -13.5%. For a better 
visuahzation the spectra are moved upwards accordingly. The field polarization is (a) parallel or 
(b) perpendicular to the rotation axis. The inset of (a) depicts a sphere and the field polarizations. 
The inset of (b) represents the cross-section variations corresponding to shape variations. 

seen in the LSPR spectrum plotted in Fig [3]d. At the same time the DF of the octopolar 
eigenmode increases with AV/V. 

C. Changes due to the variation from prolate shape 

The changes of the electrostatic resonances due to smooth variation from a prolate shape 
are shown in Fig. HJ where only the two most relevant eigenvalues and their weights are 
presented. The aspect ratio is 2. For both field polarizations, AV/V varies over a wider 
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range from -24% to 11%. The behavior of LSPR spectra with respect to shape variations 
can be seen in Fig. [51 The shape variations that are defined by volume variations can be 
vizuahzed in the left inset of Fig El Like in the previous case, the quadrupole eigenmode is 
dark, regardless of shape change. When the polarization is longitudinal, the DFs of the most 
relevant eigenmodes (dipolar and octopolar) vary smoothly over the entire range of volume 
variation considered here. The corresponding LSPRs can be noticed at longer wavelengths in 
Fig. |5l While the DF of the octopole eigenmode increases with ISV jV ^ the DF of the dipole 
eigenmode reaches its maximum at ISV jV = —11%. The LSPRs are sensitive to shape 
variations, such that at AV /V = —11%, an octopolar bump can be seen in the spectrum. 
The bump transforms into a distinct LSPR at AV /V = —23%. 

The behavior is however quite different for transverse polarization, where the cross-section 
is smaller. Dramatic changes in the weights of the relevant eigenmodes are found although 
the eigenvalues and their DFs vary smoothly. The appearance of the octopolar eigenmode 
in the LSPR spectrum occurs for a volume variation as small as ±5%. The DF of the 
dipole eigenmode is non-monotonic with a minimum at AV/V = and a maximum around 
AV/V = —11%. An interesting behavior occurs around AV/V = —10%, where the shape 
looks rather like an American football (Fig. [6^) and the weight of the octopolar eigenmode 
surpasses the weight of the dipolar eigenmode. The weight of the dipole eigenmode varies 
from 100% at AV/V = to a value of around 40% for a volume variation of -11% (Fig 
Hb). Moreover, for small positive and negative AV/V, the LSPR of the dipole eigenmode 
blueshifts as a result of the DF increase. The large reduction of the dipolar weight is 
reflected in large increase of the octopolar weight, which can be also seen in the LSPR 
spectrum plotted in Fig [5]d. At the same time the DF of the octopolar eigenmode reaches 
its maximum at AV/V = —3%. 

To understand why the octopole weight exceeds the dipole weight it is better to inspect the 
dipole and the octopole eigenfunctions of M for both shapes: the spheroidal prolate shape 
and the football-like shape obtained by a volume variation of -11%. These eigenfunctions are 
plotted in Fig. El Due to the axial symmetry, for a transverse polarization, the eigenvectors 
have the following expression u {z, ip) = Ui {z) sin when the field is along the y-axis. The 
functions Mi {z) of both the dipole and the octopole eigenmodes are plotted in Fig. [51 Thus 
the dipole eigenmode represents a large dipole (t) concentrated in the middle of the NP 
in the case of spheroidal shape (Fig. [Bb) and splits into two parallel dipoles (tt) towards 
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FIG. 4: Changes of the two most relevant DFs (left side and continuous lines) and of their weights 
(right side and dotted lines) with respect to shape variation from prolate spheroidal shape. The 
shape variation is accounted by AF/F. The first (dipolar) /second (octopolar) eigenmode is denoted 
by the square/star symbol. The field polarization is (a) parallel or (b) perpendicular to the rotation 
axis. 

the extremities for a football-like shape (Fig. [6^). These two parallel dipoles are much 
smaller than the spheroidal dipole because the extremities are much tighter than the middle 
of the NP. It is worth noticing that the dipolar character of the eigenmode is provided by 
the factor simp of the eigenfunction. At the same time, the octopole eigenmode has the 
dipole configuration 4t4 but the total dipole moment vanishes for the spheroidal shape and 
it is different from for the football shape. Also the dipolar contribution of the octopole 
eigenmode exceeds the contribution of the dipole eigenmode because the football-like NP is 
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FIG. 5: The LSPR spectra for three different AF/F: 0%, -11%, and -23%. For a better visu- 
ahzation the spectra are moved upwards accordingly. The field polarization is parallel (dotted 
lines and the LSPRs at longer wavelengths) or perpendicular (solid lines and the LSPRs at shorter 
wavelenghts) to the rotation axis. The insets depict a prolate spheroid and the field polarizations 
(right side) and the axial cross-section variations corresponding to shape variations (left side). 

much thicker in the middle than toward the extremities. 
D. Changes due to the variation from oblate shape 

The variations of the DFs and of their weights for the most relevant eigenvalues with 
respect to shape variation of an oblate spheroid and field polarization are presented in Fig. 
[71 The aspect ratio is taken to be 1/2. The axial polarization exhibits three relevant eigen- 
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FIG. 6: The z-dependent part ui(z) of the dipole (sohd hne) and octopole (dotted hne) eigen- 
functions for (a) football-Uke and (b) prolate spheroidal NPs. For transverse polarization the full 
eigenfunctions have an extra factor sin (J). The shape of these NPs are depicted in the insets. 

modes, while the transverse polarization has just two relevant eigenmodes. The variation 
of the relative volume AV /V spans between -14% and 14%. Fig. [S] shows the variation of 
LSPR spectra with respect to both shape variation and field polarization. 

The electric field parallel with the rotation axis is considered first. In this case the 
LSPRs are found at shorter wavelenghts in Fig. [HI The first relevant eigenmode represents 
also the dipole response and its corresponding DF has a non-monotonic behavior by having 
the minimum at AV/V = 4%. Thus there is a blueshift of the corresponding LSPR for 
negative AV/V (see Fig. [8]). The next most relevant is the octopole eigenmode which can 
be seen as a bump for a volume variation of -7.4% and a quite distinct new resonance at 
AV jV = —13.5%. The weight of the third relevant eigenmode does not exceed 2% over entire 
range of AV/V, therefore this eigenmode is not seen in the LSPR spectrum. The variations 
of DFs of the octopole and the third relevant eigenmode are non-monotonic. The octopole 
eigenmode has a minimum at AV/V = and a maximum around AV/V = —7.5%, while 
the third relevant eigenmode has a maximum at AV/V = 0. For transverse polarization the 
modifications of the DFs of the two most relevant eigenmodes (the dipole and the octopole) 
increase with AV/V such that for negative /positive AV/V there is a redshift /blueshift of 
the dipole LSPR. The calculations also show that the octopole eigenmode does not show up 
in the LSPR spectrum. 
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FIG. 7: Changes of the relevant DFs (left side and continuous lines) and of their weights (right 
side and dotted lines) with respect to shape variation from oblate spheroidal shape. The first 
eigenmode (dipolar) is denoted by square symbol. The other eigenmodes are denoted by star and 
diamond symbols. The field polarization is (a) parallel or (b) perpendicular to the rotation axis. 

IV. CONCLUSIONS 

In the quasistatic approximation I use a BIE method to analyze the changes of LSPRs 
produced by small and smooth variations of the shape of metallic NPs. The LSPRs are 
determined by the eigenvalues and their coupling weights that come from the operator as- 
sociated to the BIE method. First, the numerical implementation of BIE is verified against 
well known results. Thus for axially symmetric structures the method has an excellent ac- 
curacy. Also a compact formula of the NP polarizability is used to elaborate a criterion for 
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FIG. 8: The LSPR spectra for three different AV/V: 0%, -7.4%, and -13.5%. For a better 
visuahzation the spectra are moved upwards accordingly. The field polarization is parallel (full 
line and the LSPRs at shorter wavelengths) or perpendicular (dotted line and the LSPRs at longer 
wavelenghts) to the rotation axis. The inset depicts an oblate spheroid (right side) and the axial 
cross-section variations corresponding to volume variations (left side). 

the structural changes of the spectrum with respect to shape variation. Such shape varia- 
tions are inherent to NP fabrication process. Shape variations induce changes in both the 
eigenvalues and their coupling weights to the electromagnetic field. However, the coupling 
weights appear to be more affected by shape changes, such that they manifest as structural 
modifications of the LSPR spectrum. Full numerical calculations show that structural mod- 
ifications of the optical response occur for shape changes corresponding to volume variations 
greater than 12%. This quantitative result can be directly related to the criterion provided 
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by Eq. (fTT]) . where a variation of pk > 0.1 would generate structural changes in LSPR 
spectrum. A 0.1 variation of pk means a 10% volume variation via the definition of pk, 
which is inverse proportional to volume V. Therefore the NP volume variation contributes 
significantly to the structural changes of the LSPR spectrum which is agreement to other 
results regarding the averaging effect of the roughness on the LSPRs of metallic NPs^^ . The 
largest optical changes are encountered for light polarization that is parallel to the smallest 
cross-sections of NPs. This aspect is important because more care must be considered for 
the design and the fabrication of the smallest features of the NP systems. 
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